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Abstracts 

We study finite size effects of the d=3 XY model in terms of the chiral perturbation 
theory. We calculate by Monte Carlo simulations physical quantities which are, to order of 
(1/L) 2 , uniquely determined only by two low energy constants. They are the magnetization 
and the helicity modulus (or the Goldstone boson decay constant) in infinite volume. We 
also pay a special attention to the region of the validity of the two possible expansions in 
the theory. 
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1. Introduction 

One has to cope with finite size effects in study of numerical simulations on lattices. 
Recently Hasenfratz and Leutwyler have applied to this issue the chiral perturbation 
theory || which is systematic and quite alternative to the finite size scaling || [Q. The 
chiral perturbation theory is originally a low energy effective theory of QCD and it has some 
systematic expansions in momentum of pions. For a system in a finite box with volume 
V = L d , which involves Goldstone boson(s), long range properties are dominated by the 
finite mass of Goldstone bosons. By controlling its effect in a way, the chiral perturbation 
theory enables one to calculate physical quantities in a systematic series in terms of 1/L 
. As a way of controlling, an external source coupled to the field in consideration plays 
an important role. There are two manners, one is called the e-expansion, and another is 
the p-expansion, which result in the series in different powers of 1/L. Each of them has 
its own characteristic region in the space of the the external source j. The e-expansion 
characterizes the domain where j ~ L~ d , while the p-expansion corresponds to j ~ L~ 2 . 

In the present paper we consider the d = 3 dimensional 0(2) non-linear sigma model, 
or Xy-model. Though this model is simple, it is a good laboratory for studying the critical 
phenomena of the lattice scalar field theory with a global symmetry. From a realistic point 
of view as well, it is significant in connection with superfluid 4 He |5[. This model has widely 



been studied [pfl M M M laOl Ol [121 1131. It is therefore suitable to choose this model and 



to compare with those analysis in order to check whether the chiral perturbation theory 
works. The advantage of taking the model of the dimension three (d = 3) in the chiral 
perturbation theory is that finite size corrections to the partition function are uniquely 
fixed only by the two (three) low energy constants in the e- (p-) expansion not only to the 
leading order 1/L but also to the next leading 1/L 2 [JTJ] . They are the magnetization £ and 
the helicity modulus T (and an additional constant for the p-expansion) in infinite volume. 
The helicity modulus is a quantity associated with the response to the distortion of the 
direction of the spontaneous symmetry breaking |n|] . An increment of the free energy per 
volume due to the distortion is given by Af = T(Ae) 2 /2, where Ae stands for the gradient 
of the twist in one spatial direction . The equivalence between the helicity modulus and the 
Goldstone boson decay constant F is shown in ref.[|l|; T = F 2 . A numerical study of the 
finite size behavior of the system would then provide these constants and their associated 
critical indices in a high accuracy. To the d = 4 0(4) scalar model |15j |16j |L7j and the 



d = 3 classical Heisenberg model |18| the chiral perturbation theory has successfully been 
applied. 
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As stated above, each expansion has its own characteristic region, which manifests 
itself in the external source plane. It is then worthwhile to have a careful look at where 
each of them is located and whether the two overlaps or not. The authors of the paper 
15H have addressed this issue. We shall study this aspect in more detail. 

This paper is organized as follows. In section 2 we define the model and fix the 
notations. We summarize the formulae of the e- and p- expansions for the magnetization, 
the susceptibilities and the two point correlation functions. In section 3 we present the 
results of our numerical simulations based upon the cluster algorithm. Section 4 is devoted 
to conclusions and discussion. 



2. The XY model and formulae of the chiral perturbation theory in 3 dimen- 
sions 

The 0{2) linear a model coupled with a constant external source j is defined in the 
continuum euclidean space time as 

£ = l(d^(x)) 2 + ^m 2 (^(x)) 2 + A((^(x)) 2 ) 2 - jj>°(x) (2.1) 

where 4>(x) is a two component real vector {i= and 1), and m and A stand for the 
bare mass and the bare quartic coupling constant, respectively. The external source has 
only i = component. To put the model on a lattice, one rescales all the dimensionful 
quantities by a lattice constant a, and conventionally uses three dimensionless parameters; 
the quartic coupling parameter \i a t, the hopping parameter k and the external source J. 
The lattice action is 

S = ~ k E + E titi + A ^ E (E titi - 1) 2 - j E (2.2) 

n,/i,i n,i n i n 

where (p l n is the dimensionless field <p % n = </> l (na)a d / 2_1 sitting at site n and \i is a 

unit vector of the // direction. The lattice parameters are related to those of continuum 
in such a way as (ma) 2 = (2 — 4A^ at )/K — 2d, Xa 4 ~ d = Xi a t/^ 2 , ja 1+d / 2 = J/y/n. For 
Xi a t = oo, the radial mode of tp n is frozen to unity, and the action reduces to that of the 
0(2) non-linear sigma model or the XY model. 

Our aim in the present paper is to compute, in the d = 3 XY model, two low energy 
constants in infinite volume, which are the magnetization S 

hm hm (0%, = £ (2.3) 

j— >0 V -^oo 
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and the helicity modulus (or the Goldstone boson decay constant) F appearing in the 
following equation 



J dx e ipx (A^x)A u (0))^o,v=oo = ^ v% ^ + • • • (2.4) 

where = (jPd^cf) 1 — (j^d^cjP . The quantity {(j>°)j,v is the expectation value of <p° in a 
three dimensional box with finite volume V in the presence of the external source j. 

In d = 3, all the formulae of the chiral perturbation theory are basically fixed only 
by these two low energy constants to 0(1/L 2 ). It is unlike the four dimensional models 
where a few additional constants are necessary [II. In the following, we summarize the 
formulae of the chiral perturbation theory necessary to our analysis of finite size effects 
on the magnetization, the susceptibility and the two point functions. Detail about the 
formulae is found in ref. for the e-expansion and in ref. |19|] for the p-expansion. 

In order for the chiral perturbation theory to be effective, the Goldstone modes should 
be important at long distances. That is to say, finite size effects from the massive compo- 
nent with a mass m a must be negligible compared to the Goldstone boson mass m n , i.e., 
m cr /m 7T 3> 1. Assuming the dominance of the Goldstone modes two manners of expan- 
sions are possible depending upon what quantities are fixed in the expansions One is 
called e-expansion, in which the total magnetic energy uq = T,jL d is fixed to O(L ), and 
another is p-expansion where v = T,jL 2 /F 2 is fixed to O(L ). The latter gives j ~ L , 
while the former does j ~ 1/L d . By using the Goldstone mass m^, which is given by 

= Ej/F 2 , or the corresponding correlation length £ w = l/m f , the domain is char- 
acterized by m^L < 1 (£ n > L) for the e-expansion and by m^L > 1 (£ w < L) for the 
p-expansion. Since m^L = ^/YPjLjF < 1 yields j < F 2 /T,L 2 = const, for a fixed k and L, 
one expects in the j space that the e- (p-) expansion holds in the smaller (larger) j domain. 
As one moves k, each of the regions shifts. In the vicinity of the critical point k c , E and 
F behave as E ~ (k — k c )"( 1+v ^ 2 and F ~ (k — k c ) u I 2 , respectively. One hence expects 
that m^L = ^/ZjL/F w 1 leads to j ~ F 2 /E ~ (k - K c ) iy(1_7 ' ) for a fixed L. Therefore 
the domain for the e-expansion shrinks as k approaches k c from above, since the index 
77, anomalous dimension of the field, is smaller than unity for the model in consideration. 
It is one of our motivations to locate the region of each expansion in the j space. It is 
also interesting to study to what extent each region is extended beyond the point where 
m n L w 1. This is another motivation of our work. In the present paper we shall look into 
these aspects in detail. 
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All the formulae in this section are of order 0(L~ 2 ). Hereafter we call the direction 
of the external source (i = 0) the longitudinal direction and another direction (i=l) the 
transverse one. They are denoted as || and _L, respectively. We work only with finite 
volume and finite j, and drop, for simplicity, the suffices j and V from ((jp)jy. As an 
expansion parameter we use a(= 1/F 2 L). 

2.1. e-expansion 

The magnetization in the direction of an external source is 

(0°) = Eu[p ir] + 2p 2 a% (2.5) 

where u = p\T,jV, and p\ and p 2 are quantities depending on the shape of the box; 

p 1 = l + l-(3 1 a+l((3 2 1 - 2(3 2 )a 2 

1 2 8 (2.6) 

P2 = 

For a three dimensional symmetric box, f3\ = 0.225785 and P 2 = 0.010608. The quantity 
i] in fl2"T5| ) is given by the modified Bessel's functions I n (x) as 

u I {u) 

Two point correlation functions are defined as 

^'(t) = -^^(0>,t)^(O,O)) (2.8) 

n 

where the summation is taken over in spatial 2- volume at fixed "time" t (0 < t < L). We 
consider two correlation functions; the longitudinal and transverse ones; 

Gl(r)=G 00 (t), Gl(r) = G 11 (t), (2.9) 

where r = t/L and < r < 1. In the e-expansion each correlation function takes the 
following form 

Gjj(T) = a,| +&||MT) + C||MT) + d||MT), (2-10) 
G1(t) = a ± + 6x/ii(r) + c ± h 2 (r) + d±h 3 (r), (2.11) 
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where coefficients oy , 6y , cy and d\\ and their transverse counterparts are expanded in powers 
of a. The r dependence enters through the foilowing kinematic functions 

*.W = i((r-i)»-i) 

1 '.2„ ^2 U 



MT ) = _^( 1 _ T )^__ j (212) 

■ n , , ^' cosh(g n (r- 1/2)) 

The summation on the r.h.s. of ^13 (t) is taken over all integers (m, n 2 ) except for n = (0, 0), 
and stands for g n = 27TA/n 2 + n\. The coefficients a, b, c and ei for the longitudinal 
correlation function in ( |2.1U| ) are given by 

ay = E 2 {p\(l - V ) + 2p 2 (l + 2^ 2 77)a 2 } 
6|l = E 2 p?f7« 

(2.13) 

C|| = E 2 (l-2?7)a 2 
d|| = E 2 (l -??)a 2 /2, 
while the transverse ones in ( [2.11| ) are given by 

a± = T? {p 2 r ] + 2p 2 a 2 } 

b ± = X 2 pj(l- V )a 

(2.14) 

c± = E 2 (-l + (2 + w 2 )?7) a 2 
d± = E 2 ?7a 2 /2. 

Note that a, 6, c and d are of O(l), O(a), 0(a 2 ) and 0(a 2 ), respectively, for both the cases. 
We use these correlation functions to extract the low energy constants. In addition to 
these, we find that the scalar product correlation function is also practically of much use; 

^(r) = i^(n,Lr).#,0)) 

n 

= G\\(t) + GI(t) (2 - 15) 

= a s + b s hi(r) + c s h 2 (r) + d s h 3 {r), 
where the coefficients a s etc. are given by 

a s = E 2 {pj + 4p 2 (l + u 2 7 1 )a 2 } 
b s = Yj 2 p\a 

r 2 2 2 V^-16) 

c s = h u rja 
d s = E 2 a 2 /2. 



5 



The longitudinal susceptibility is calculated from it's correlation function 



xii = E «<An)0°(o)> - <</>°> 2 ) 



(2.17) 



and it takes the form 



XII = E 2 V [p 2 {(l-r})-u 2 V 2 } + 2p 2 a 2 ] 



(2.18) 



2.2. p-expansion 

In the p-expansion to 0(L~ 2 ) a low energy constant ko gets involved in addition to E 
and F for the magnetization and the susceptibilities. 
The magnetization reads 

(0°) =£ {l + ^«C 1/2 + k a 2 C -\a[l+ ^aC 1/2 ]^i(C P ) 



+ 5 aVC)MC)-fo(C)] 



(2.19) 



where £ = (m w L) ( = T,j/F ) , and £ p is 

Cp = C(i + ^C 1/2 )- 

The function g n (C)( n = 0, 1, 2, ...) are defined by 

A3/2 00 o 

^o(C) = -lnC-^-E^ 



4n 



d 



n=0 



9n+l(0 = --jf 9n(C)- 



The correlation functions defined in (|2.9|) and (|2.15|) take the forms 



2 1. 



Gf(r) = (0°> +-E 2 a%(r) 
^(r) = Z„a/n(r), 



(2.20) 



(2.21) 
(2.22) 



(2.23) 
(2.24) 



where Z„ is the wave function renormalization constant in a finite box with the form 

(2.25) 



Z v — X' 



l + -3-aC 1/2 -^i(C)« 

47T 
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The kinematic functions hi(r) and h 2 {r) involved in the p-expansion are 

hi(r) = h(rXv) 

h 2 (r) = J2h 2 (rX + ^ 2 \n\ 2 ) 



where h on the r.h.s. is given by 



h(r, C) 



1 cosh (ygr - 1/2)) 
2~vT sinh 



and is 



1 - ^i(0« 



(2.26) 



The transverse correlation function has a normal form of massive scalar propagator to 
O(a), while the connected part of the longitudinal one ( p.23|) is 0(a 2 ). We also consider 
the scalar product correlation function 



As in the e-expansion, G p s (t) is of practical use. 

The susceptibility in ( |2.17|) in terms of p-expansion reads 



(2.27) 



ZVal—a( 1/2 + k a< 
{ 16n 



a 

2 



-92(( P ) - ^ a C 1/2 (gi(C P ) + 2(1 + C)<7 2 (Cp)) 



+ - A a 2 [-2g l92 + C(g 2 2 + gig 3 )] i . 



(2.28) 



3. Numerical simulations 

We employ, to generate configurations, the Wolff's single cluster algorithm |2(| and its 



modification in the case of presence of an external source|T8[. The lattice shape which we 
take is a symmetric box with a volume V = L 3 , and the size of the box ranges from L = 32 
to L = 64. The d = 3 dimensional XY model is known to have a second order phase 
transition at k c = 0.45420(2)0. We calculate at k larger than k c , and choose two k values 
K—0A62 and 0.47 for main calculations of the low energy constants, and several more for 
obtaining the critical indices. The magnitudes of the external source changes from J = 0.0 
to J = 5.0 x 10 -4 . In terms of continuum external source, which depends on k, j ranges 
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from to 7.29 xlO" 4 for, say, k = 0.47. We make measurements at each 5 updatings with 
typical number of configurations for measurements ranging from about 20,000 to 40,000. 
The statistical errors are estimated by using the blocking and the bootstrap methods. We 
use the program system SALS |yj for the least square fitting of the data to the formulae. 



We extracted £ and F from three types of correlation functions in each of the e- and 
p- expansions. We looked at the behaviors of each of the extracted values by varying J. 
We also observed the magnetization and the susceptibility, and compared the results with 
the formulae of the chiral perturbation theory. Main results are in order. 

(1) We found, as expected from the arguments in the Sec. 2, that the e-expansion provides 
good results in the small J region, while the p-expansion does so in the larger J region. 
We identified where the region of the validity for each of the two manners of expansions 
is located. 

(2) The scalar product correlation functions provide stable results for the two constants. 

(3) Upon varying the size of the lattice, which is sensitive to the region of the validity, 
we observed a significant change of the behavior for the p-expansion. It is consistent 
with the arguments concerning m n and m a . 

(4) The magnetization shows a clear crossover behavior between the two expansions as a 
function of J. The longitudinal susceptibility turns out too poor to address the region 
of the validity. 

(5) We calculated the critical indices associated with the two constants. The results are 
consistent with other references. 

In the following we shall discuss the results in detail. 

3.1. Correlation functions 

In order to extract and F(k) from the correlation functions, we fit to the for- 

mulae (FOCI) , O and (ETI5D for the e-expansion and to (|2~23l) , ($T2§ and flT27D for 



the p-expansion for each j at fixed k. We use hereafter the notation ^(e,p),(\\,±,s) f° r the 
value of E(k) obtained from each type of correlation functions in each expansions. Similar 
notations are used also for F. All the extracted values of £ and F are listed in Table 
I - V. In what follows we discuss the stability of the extracted results and the region of 
the validity of each of the expansions by comparing the values in Tables. If the chiral 
perturbation theory works, the extracted values should be independent of J, and all of the 
values £( c ,p),(||,_L, 8 ) 0P( c ,p),(||,j_,.)) should agree. 



Results of L = 32 lattice (k = 0.47) 

In Figs. 1 and 2 we show the behaviors of the obtained £(«) and v.s. the 

magnitude of the external source J at L = 32 for k=0.47. In Fig. 1 we see that the 
e-expansion provides the consistent result for E within errors for two types of correlation 
functions, i.e., longitudinal one Gjj(r) and transverse one Gj_(t), in the whole J range 
in consideration. The p-expansion, on the other hand, is out of validity in the smaller J 
region, J < 2.0 x 10~ 4 (Fig. 1). At larger J values, the transverse one E Pi _l agrees with 
£ of the e-expansion. We hence observe that an overlap of the two expansions takes place 
at J > 2.0 x 10" 4 . 

The result of the longitudinal p-expansion is not included in the figure. Fittings 
according to this type of expansion turn out invalid in the whole J range. This may be 
due to the fact that r dependence appears only at 0(a 2 ) unlike the all other correlation 
functions of 0(a), and that the data would be too noisy to extract the E. This feature 
applies also to F and to different volumes. We then drop the result of the longitudinal 
p-expansion throughout the paper. 

Fig. 2 shows F v.s. J. Both of F e ^ and F £) j_ agree within errors in the small J region 
(J < 2.0 x 10~ 4 ), while F e || and F C: ± split from each other for larger J (J > 2.0 x 10~ 4 ). 
This indicates that the boundary of the two domains where the e-expansion is valid and 
invalid is located at J « 2.0 x 10~ 4 . Fig. 2 also shows the behaviors of F for the 
p-expansion. Similarly to E, the p-expansion does not give a reasonable result in the 
smaller J region, but agree with F of e-expansion at J > 2.0 ~ 3.0 x 10~ 4 . 

The lowest value of the overlapping region for E and F is consistent with what one 
naively expects from the argument concerning the Goldstone mass (Sec. 2); the two 
expansions coincide at a point J which is determined by the relation M n L = ^/EjL/F ps 1. 
Taking the values £(0.47) = 0.2581(1), F(0.47) = 0.286(1) and L = 32 for k = 0.47, it 
yields 2.1 x 10 -4 , which is in good agreement with the observed value. Note that the errors 
of F are about five times larger than those of E at whole k's and J's in consideration. 

Results of L = 32 lattice (k = 0.462) 

Similar behaviors are seen for the data of k = 0.462 and L = 32 (Fig. 3). However, E e y 
and E £) _|_ agree only at small J values (J < 2.0 x 10 -4 ) compared to k = 0.47 (Fig. 1). 
This is in fact what is expected from the condition giving the region of the expansion 
M^L < 1, since the region of the e-expansion shrinks according to j < (k — k c ) i/< ^ 1_?7 - ) as 
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K approaches k c (k c = 0.4546) as stated in Sec. 2. We see that this behavior becomes even 
clearer for F (see Table II). 

Fig. 3 also shows the behaviors of the p-expansion. The value E p ,_i_ seems reasonable 
only in the limited intermediate region around J ~ 2.0 x 10 -4 . In larger J region the 
p-expansion becomes invalid unlike the case of k = 0.47. The reason for that comes 
from the mass m a of the massive particle; as stated in Sec. 2, m a must be large enough 
rricr/m-n- ^> 1. In the case of k = 0.47 we obtained m a ~ 0.41, which gives the ratio ma/rn^ 
ranging from 20.5 to 10.1 as J changes from 1.0 x 10 -4 to 5.0 x 10 -4 . m a is then heavy 
enough for the chiral perturbation to apply at k = 0.47. 

For k = 0.462, on the other hand, we obtained smaller value m a = 0.3, which yields 
rrifj/m^ ranging from 12.5 to 5.5 in the same J range. The authors of the paper [^5j took 



rria/m^ > 10 as a criterion. If we employ the same value, it provides an upper bound of J, 
J < 1.5 x 10 -4 , on the region of the validity of the chiral perturbation theory. Combining 
the condition m^L > 1, which gives a lower bound of J, J > 1.66 x 10 -4 , on the region of 
the p-expansion, one observes that there is no region for the validity of the p-expansion. If 
m a L > 7.5 is adopted instead of 10, however, the upper bound shifts to J = 2.7 x 10~ 4 . 
Although these bounds should not be taken strictly, it indicates that the p-expansion is 
expected to apply to no or only the narrow J region around J = 1.6 ~ 2.5 x 10 -4 . In 
other words, the "window" for the p-expansion is almost closed or open only slightly. The 
results of our fittings are in good agreement with this inference (Fig. 3). 

Scalar product correlation function 

Apart from the longitudinal and transverse correlation functions, we made also use 
of the scalar product correlation functions ( |2.15| ) and ( |2.27| ). During our fitting procedure 
we came to notice that the scalar product correlation functions are good correlators which 
provide stable fitting results. In particular, the e-expansion works very well. 

In Fig. 4 we show the results of £ for both of the e- and p-expansions at L = 32 for 
k = 0.47 and 0.462. For both k values E e jS looks almost independent of J in the whole 
J range, and therefore it provides reliable result. This is the reason why we indicated the 
value E £jS in Figs. 1 and 3 by the arrow as a reference. In the region of validity of the 
e-expansion at each k, E e s is consistent with E e m and £ ej j_. Compared to the other two, 
E e s is much stable. 

The value E PiS , on the other hand, vary as J moves. For k = 0.47 it monotonically 
approaches E e>s as J increases, and gives consistent result at J > 3.0 x 10~ 4 . For k = 0.462, 



10 



on the other hand, E PjS never agrees with E £jS even in the larger J region. The scalar 
product correlation functions supports more clearly, than do the longitudinal and the 
transverse correlation functions, the inference based upon m a and about the location 
of the boundary between the two expansions. Similar behaviors are seen also for F as 
shown in Fig. 5. 

Volume dependence 

Let us turn to the volume dependence. The extracted values of E and F for L — 48 
are listed in Tables III and IV . The average values for L = 48 are consistent with those 
for L = 32 in the regions of the validity of each of the expansions. 

When one increases the volume, J value corresponding to the condition m^L rs 1 
decreases as J oc 1/L 2 , while the ratio m^/m^ is independent of L. It is therefore expected 
that the "window" for the p-expansion becomes wider as L increases. It is then interesting 
if one actually sees this behavior by simulations. For L = 32 at k = 0.462, we have seen 
that the "window" is almost closed. Fig. 6 shows E £jS and E PjS for L = 48 at k = 0.462 
obtained from the scalar correlation functions. The value E P)S agrees with E 6)S only around 
J = 2.0 x 10 -4 , which is regarded as an evidence of slightly opening "window". It is in 
good contrast with the behavior of E(k = 0.462) in Fig. 4. Similar behavior is found also 
for F as shown in Figs. 5 and 7. 

In the absence of J we calculated E and F at L = 64, and compared them with those 
of L = 32 and 48. As seen in Table V the results of E £jS and F BiS are consistent within 
errors. 

3.2. Magnetization and susceptibility 

We calculate the magnetization for the e-expansion ([?75|) for the p-expansion ( |2.19| ) . 
Rather than fitting ( |2.5[ ) and ( |2.18| ) to the data to extract E and F, we put E and F into 



eqs. ( |2.5| ) and ( |2.19| ) and compare them with those of direct measurements of (<p°). Here 
we use the values of E and F obtained from the correlation functions at each k, k = 0.462 
and k = 0.47, in the previous subsection. In the p-expansion, an additional low energy 
constant ko is involved as mentioned before. To determine ho, we match the curve ( |2.19| ) 
to the data at some J. 

Fig. 8 shows the results for k = 0.462 and 0.47 at L = 32. The constant ko is fixed at 
J = 2.0 x 10" 4 for k = 0.462 and at J = 3.0 x 10" 4 for k = 0.47. Solid line in the figure 
shows the curve (|2.5| ) and dotted line is ( |2.19| ). We see a good agreement and observe a 
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clear crossover between the e- and the p-expansions. The location of the crossover between 
both the expansions is consistent with the prediction m^L pa 1. 

The results of L = 48 are shown in Table VI. The values obtained from the p-expansion 
formula ( |Tg ) are (<p°) = 0.201 at J = 0.5 x 10" 4 , 0.2934 at J = 2.0 x 10" 4 and 0.3083 
at J = 3.0 x 10" 4 for k = 0.462. The two of them (J = 0.5 x 10" 4 and 2.0 x 10" 4 ) 
are consistent with the numerical result within errors. Another one ( J = 3.0 x 10 -4 ) 
seems to deviate slightly from the data. For k = 0.47, the values of (if ) read 0.279 at 
J = 0.5 x 10" 4 , 0.3357 at J = 1.0 x 10" 4 and 0.3776 at J = 3.0 x 10" 4 , respectively. In 
the similar manner to « = 0.462 only the two points (J = 0.5 x 10 -4 and 1.0 x 10 -4 ) 
out of three are in good agreement with the data within errors. In any case, theoretical 
predictions are not inconsistent with the measurements. 

Let us turn to the longitudinal susceptibility x\\- Fig. 9 shows the J dependence of 
X|| for L = 32 at k = 0.462 and 0.47. The solid (dotted) lines correspond to (^T8| )( (l2~28D ). 
For each of the fixed k, x\\ shows a crossover at J ~ 2.0 x 10 -4 if the e- (p-) expansion 
is valid in the smaller (larger) J region. In the actual measurements we are not able to 
observe such a clear crossover. 

As to the relation between the lines for the two k's we can read off qualitative difference 
as follows. For smaller J (J < 2.0 x 10 -4 ), the values of X\\{ K = 0.462) are lower than 
X\\( K = 0.47), while in larger J region X\\{ K = 0.462) becomes higher than X\\( K = 0.47). 
This property is successfully seen in the direct measurements. Quantitatively, however, 
there are some inconsistencies between the theoretical predictions and the numerical data. 

The J dependence of X|| in the case of L — 48 is shown in Table VII. The situation is 
the same as in the case of L = 32. 

3.3. Critical index 

The critical indices f3 and v are defined by 



in 3 dimensions 0. To extract f3 and v we make use of the values of E and F in infinite 
volume obtained from the previous analysis. We fit to (|3.1| ) the data at five k points 
(k = 0.46, 0.462, 0.464, 0.466, 0.47). The results are 



S(k) ~ {K-K c f 

F(«) ~ (k - K c y' 2 



(3-1) 



k c = 0.4546(5) 



(3.2) 
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for the critical coupling 

13 = 0.321(17) 
v = 0.66(6) 

for the indices. They are consistent with other references! 



(3.3) 



™rap]ra. 



4. Conclusions and discussion 

We applied the chiral perturbation theory a la Hasenfratz and Leutwyler to the d = 3 
XY model in order to calculate the two low energy constants. They are the magnetization 
E and the helicity modulus F (or Goldstone boson coupling) in infinite volume. In the 
theory, two manners of the expansions are involved. One is the e-expansion, which is 
valid in the region where m n L < 1, and another is the p-expansion, where m^L > 1. On 
L = 32, 48 and 64 lattices, we fitted the formulae of the correlation functions to the Monte 
Carlo data. All the values of E and F extracted in each of the regions of the validity are 
consistent and volume independent within errors. 

We are also particularly concerned with the region of the validity of the two expansions. 
As far as the two k values are concerned, the lower boundary of the p-expansion is basically 
located at the region where m^L pa 1 holds. On the other hand, the e-expansion stretches, 
for some cases, beyond the point M^L pa 1 to some extent than expected. The similar 
behavior was also observed in ref.[15|]. Apart from the condition for m ff L, the mass m a of 



the massive particle puts a constraint on the validity of the chiral perturbation theory. As 
expected from these conditions, we have observed the significant difference of the behaviors 
of E and F v.s. the external source J for the two different k values. 

We found that the scalar product correlation functions are better correlators than the 
longitudinal and transverse correlation functions. Particularly the e-expansion provides 
quite stable estimations of E and F. By use of it the overlapping of the validity of both 
the expansions is clearly seen. The reasons for the stability are in order. The coefficient 
b s = E 2 p\a in ( |2.16|) is independent of j, unlike the longitudinal bu in ( |2.13|) and the 



transverse counterparts b± in ( |2.14j ), in which j dependence appears through u = piT,jV. 
To 0(a), therefore, no response of G e s to the variation of j appears. It is then expected 
that E e s and F e>s are extracted independently of j. This is one of the reasons for the 
stability *. However for J > 3.0 x 10 -4 , E and F deviate from their stable values, i.e., in 



The fact that b s is independent of j to 0(a) applies only to N = 2 for the O(N) model. 
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this J region the contribution from 0(a 2 ) would be responsible for the stability. In order 
to have a look at its effect, we compared the two ways of fitting of G% using the formulae to 
0(a) and to 0(a 2 ). We found that in the small j region (J < 3.0 x 10 -4 ) the results of the 
two fittings are in agreement with each other within errors, while in the larger region the 
significant discrepancy appears. (The relative difference of the two fits is approximately 
4% for F(k = 0.47) at J = 5.0 x 10~ 4 ) One then sees that the correction of 0(a 2 ) plays a 
significant role to the stability in the larger J region. 

The magnetization ((p°) and the susceptibility x\\ are fitted by use of the values of £ 
and F obtained by the correlation functions. The consistency was checked. The magne- 
tization is well fitted by the e-expansion in the smaller J region and by the p-expansion 
in the larger one. A crossover between both the expansions in the J space was observed, 
and its location is consistent with the prediction m^L w 1. The susceptibility turns out 
too poor in precision to address the issue about the validity of the expansions. 

We calculated the critical indices associated with S and F, and obtained the values 
consistent with other references. However we have not reached the precision as high as 
that of ref.'s || and ||, which is based upon the phenomenological renormalization group. 
This may be due to the fact that the chiral perturbation theory becomes hard to apply as 
it is. Because, as one approaches the critical point, the mass m a of the massive particle 
gets smaller and the ratio rria/m^ accordingly becomes small. In addition, the expansion 
parameter a = 1/F 2 L becomes large. 

A few words for the values of the constant F. We compared the value of T calculated 
from F with the one in the available references. Our values** T = F 2 = 0.0515(14) at 
K = 0.462 and 0.042(2) at k = 0.460 in infinite volume come slightly below the curves in 
the figures of refs. || and [||, in which the lattice size L is at most 48. It is reasonable 
since in their data the values monotonically decrease as L grows. 

** Note that F 2 / k corresponds to the helicity modulus in their literatures. 
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Table Captions 

Table I E and F for k = 0.47 and L = 32. The symbol *** indicates the point where the 
p-expansion is beyond the applicability. The symbol ### shows that we failed to fit 
the formulae. The same symbols are used in the other tables. 
Table II E and F for k = 0.462 and L = 32. 
Table III E and F for k = 0.47 and L = 48. 
Table IV E and F for k = 0.462 and L = 48. 
Table V E and F for L = 64 and J = 0.0. 

Table VI ((p°) v.s. J for L = 32 and L = 48. By the symbol %%% it is meant that we have no 

data available. 
Table VII x\\ v.s. J for L = 32 and L = 48. 
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Figure Captions 

Fig. 1 E v.s. J for k = 0.47 and L = 32. Circle indicates E e ||, square is E €; j_, and triangle 
corresponds to E p> x- For small J there is a significant difference between the e- 
and p-expansions. The arrow shows the mean value of E £jS , as a reference, which is 
calculated from E £jS at seven J points in Table I. The value E £)S is regarded as the 
best estimate of E. Detail about E £)S is found in the latter part in this section. 

Fig. 2 F v.s. J for k = 0.47 and L = 32. Circle indicates F e u, square is F e> ±, and triangle 
corresponds to F P) ±. The location of the arrow shows the mean value of F e>3 , as a 
reference, which is estimated in the same manner as in Fig. 1 from Table I. It is also 
the best estimate of F. 

Fig. 3 E v.s. J for k = 0.462 and L = 32. Circle indicates E e n, square is E £) j_, and triangle 
corresponds to E P) j_. For J > 2.0 x 10 -4 the e-expansion appears out of validity. The 
arrow shows the mean value of E £ S . 

Fig. 4 E v.s. J for L = 32. All of the values E are determined by the scalar product 
correlation functions. E £ S is almost independent of J for both the k values. For 
k = 0.47 E P)S (square) agrees with E £iS (circle) at J > 3.0 x 10 -4 , while for k = 0.462 
there is a significant difference between E £ S (filled circle) and E p s (filled square). The 
statistical errors of E £ S lie within the symbols. 

Fig. 5 F v.s. J for L = 32. All of the values F are determined by the scalar product 
correlation functions. Similar behavior to E in Fig. 4 is observed. The same symbols 
as those in Fig. 4 are used. 

Fig. 6 E determined by the scalar product correlation functions at L = 48 and for k = 0.462. 
E PjS (square) agrees with E £iS (circle) at J = 2.0 x 10~ 4 . As seen in Table IV, the 
p-expansion does not provide reasonable value for E p s at J = 0.5 x 10~ 4 , i.e, out of 
validity. The arrow shows the mean value of E £ S . 

Fig. 7 F v.s. J at L = 48 and for k = 0.462. F P;S (square) agrees with F £jS (circle) at 
J = 2.0 x 10 -4 . The arrow indicates the mean value of F e>a 

Fig. 8 Magnetization (ip°) v.s. J for L = 32. Normalization of the data is done at J = 
3.0 x 10" 4 for k = 0.47 and at J = 2.0 x 10" 4 for k = 0.462. The solid lines show 
(F2~5|) (e-expansion) , and the dotted lines correspond to (|2. 19|) (p-expansion) . A clear 
crossover between the two expansions is observed for both k = 0.47 (square) and 
k = 0.462 (circle). 

Fig. 9 Longitudinal susceptibility x\\ v - s - J for L = 32. The solid lines show (|2.18|) (e- 
expansion), and the dotted lines correspond to ( |2.28| ) (p-expansion) . 
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